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We study the transport properties of disordered two-dimensional electron systems with a 
perfectly conducting channel. We introduce an asymmetric Chalker-Coddington network model 
and numerically investigate the point-contact conductance. We find that the behavior of the 
conductance in this model is completely different from that in the symmetric model. Even in 
the Umit of a large distance between the contacts, we find a broad distribution of conductance 
and a non-trivial power law dependence of the averaged conductance on the system width. Our 
results are applicable to systems such as zigzag graphene nano-ribbons where the numbers of 
left-going and right-going channels are different. 
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The critical behavior of the transport properties of 
two-dimensional electron systems under quantum Hall^^ 
conditions has been investigated in various models. The 
Chalker-Coddington (CC) network model^^ (Fig. 1) is 
especially suited to the calculation of transport proper- 
ties^^ because current amplitudes are calculated directly. 
This is in contrast to the tight binding model where the 
wave functions must be calculated first. 

The CC model consists of links corresponding to 
equipotential lines and nodes describing the scattering 
at saddle points of the random potential. Assuming the 
amplitudes of incoming and outgoing currents for a node 
to be ci,C4 and C2,C3, respectively (see Fig. 2), scatter- 
ing at a node is described by a 2 x 2 unitary scattering 
matrix s, 











(1) 



(2) 



The effect of disorder is included in the phases (/)^, which 
are independently and uniformly distributed between 
and 27r. The scattering probability p controls whether 
the system is in the insulating regime {0 < p < 0.5), at 
criticality {p = 0.5), or in the quantum Hall insulating 
regime (0.5 < p < 1). 

1.1 Asymmetric Chalker-Coddington network model 

In tight binding models, the numbers of the right- 
going and left-going channels are always the same. Under 
certain conditions, however, some of the channels decou- 
ple. One example is a graphene sheet with zigzag edges, 
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(a) insulator 
(P^O) 



(b) criticality 
(p = 0.5) 



(c) quantum Hall 
insulator (p^l) 





Fig. 1. A schematic of the current flow in the CC model. Current 
flows along the links: arrows indicate the direction of flow. In 
the insulator (a), all current paths are closed and all states are 
localized, while at criticality (b), all states are delocalized. In 
the quantum Hall insulator (c), bulk states are localized, but the 
edge states carry current. 




Fig. 2. Scattering at a node. The current amplitude on link i is c^. 
The incoming currents are scattered to the left with probability 
p and to the right with probability (1 — ^>). 

where there are say, left-going and / + 1 right-going 
channels near /ca = 27r/3, and / + 1 left-going and / right- 
going ones near ka = — 27r/3, where k is the wave num- 
ber and a the lattice constant. For long ranged scatter- 
ers, states near 27r/3 and — 27r/3 do not mix, and hence 
the numbers of right-going and left-going channels be- 
come, in effect, asymmetric. This asymmetric situation 
has been studied numerically for quantum railroads^^ 
and analytically^' on the basis of the DMPK equa- 
tion.^' Here we realize such an asymmetric situation 
in the CC model^o) (Fig. 3). 
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(a) asymmetric 
insulator (p^O) 



(b) criticality (c) asymmetric 
{p — 0.5) insulator (p^l) 




Fig. 3. A schematic of the current flow in the asymmetric CC 
model. In the insulating phases (a) and (c), there is a conducting 
channel at one of the edges. 




Table I. The two-terminal conductance in quasi-one dimensional 
symmetric (n^ — = 0, see Fig. 1) and asymmetric (n^ — 
nout ^ Pig 3) models. 





p^O 


p = 0.5 




Symmetric 


(^qlD ^ 






Asymmetric 


Gra = l 


GrR = l 


GrR = l 



For asymmetric systems with two-terminal geome- 
try, terminals at the ends of the system are also asym- 
metric in the numbers of incoming and outgoing channels 
and the two-terminal conductances measured with cur- 
rent flowing left to right Gl^r, and right to left Gr^l 
are related by 

Gl^r = Gr^l + - nr ), (3) 

where G is measured in units of e^//i. Here, L and R refer 
to the left and right terminals, and nj^ and n^^^ are the 
number of incoming and outgoing channels, respectively, 
in the left terminal. It follows from current conservation 
that 

nL" + nk"=nr +nr- (4) 
Using this equation, we can rewrite eq. (3) as 

Gr^l = Gl^r + (ng - n^). (5) 

If we suppose that nj^ > n^^^ it follows that 

Gl^r > nL" - (6) 

Thus we expect Gl^r to be finite even in the limit of 
infinite length (see Table I). The analysis of the transmis- 
sion eigenvalues shows that the system has nj^ — ^l^* per- 
fectly conducting channels.^' ^' However, the formula 
(6) makes it appear that this property is a consequence of 
the asymmetry of the terminals rather than the sample. 

In this paper, we calculate the point-contact con- 
ductance Gpc of an asymmetric CC network model. The 
point-contact conductance is the conductance measured 
between two interior probes. Just like the probes of 
a scanning tunneling microscope, the probes make con- 
tact with the sample at a point. The probes work as 
symmetric terminals (nj^ = ^l^* = = ^r^* = 1) and 
Gpc varies between and 1 in units of e^//i. In the next 



Fig. 4. A schematic of the asymmetric CC model with point- 
contacts 

section, we explain how to calculate the point-contact 
conductance. In § 3, we show that the asymmetry of the 
network is reflected in a broad distribution of Gpc with 
finite averaged values in the long distance limit. In the 
final section, we summarize and conclude. 

2. Method 

We denote the numbers of links in the x and y di- 
rections by Lx and Ly^ respectively. We impose periodic 
boundary condition (PBC) in the x direction and fixed 
boundary conditions in the y direction. This corresponds 
to a ring geometry. For PBC in the x direction, must 
be even. We regard the links at x = + 1 as the ones at 
X = 1. In the standard CC model, Ly is even and the sys- 
tem is symmetric. Here we set Ly odd so that the system 
is asymmetric. The state of the network is specified by 
the complex current amplitudes Ci on the Lx x Ly = N 
links. 

2.1 Point-contact conductance 

To introduce point-contacts into the network, 
we cut link L at (xi^yi) and link R at (x2,^2)- We then 
define incoming current amplitudes cf^ , cg^ and outgoing 
current amplitudes Cj^^^c^^ on the corresponding hnks 
(Fig. 4). The current amplitudes satisfy the equation 















^out 


= S 


Jn 


^out 




Jn 






\cn/ 



where S is the N x N scattering matrix consist- 
ing of 2 X 2 scattering matrices s at each node. 
For given (cl,cJ^), the remaining current amplitudes 
(ci, C2, • • • , Cl^*, • • • , Cp^*, • • • , Cat) are uniquely deter- 
mined by the following set of simultaneous linear equa- 
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tion with N unknowns 







/ci\ 







^out 


-s 
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Jn 
^L 


^out 
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(8) 



As a consequence of the structure of these equations, 
there is a hnear relationship between the incoming and 
outgoing current amphtudes 



L _ 



(9) 

^R / \^ ' / \^R/' 

The most straightforward way to calculate the transmis- 
sion coefficient is to set 



so that 



t 



(10) 



(11) 



(12) 



The point-contact conductance Gpc is given by 

in units of e^/h. 
3. Results 

3.1 Distribution of point- contact conductance 

The point-contact conductance depends on the posi- 
tions {xi^yi) and (^2,^2) of contacts in addition to the 
parameters of the network L^, L^, and 



Gpc = Gpc{xi,yi,X2,y2,LxjLy,p). 



(13) 



This is a sample dependent quantity. If we average over 
disorder, translational symmetry is recovered, and the 
averaged conductance (Gpc) depends only on the dis- 
tance \xi — X2I = Xd for fixed yi and y2 (see Fig. 5). 
Taking yi=y2 = yp, (Gpc) is a function of Xd, 2/p, L^, 
Ly, and p 



(Gpc) = F{xd,yp,L^,Ly,p). 



(14) 



For convenience, we consider only two values of 7/p, cor- 
responding to edge conductance Gp^; yp = 1 and bulk 

conductance Gp^; yp = ^^2^ • 




Fig. 5. The geometry of our modeL Coordinate x is measured 
along the ring and y across the ring. The circumference of the 
ring is Lx and the width is Ly . The distance between the contacts 
is (< Lx/2). 



Table IL The point-contact conductances in the insulating limits 
(see Figs. 1 and 3). Gpc is unity only when the contacts are 
directly attached to the edge states. Otherwise, Gpc = 0. 



P - 



Symmetric 
Asymmetric 







VpLy 



In the insulating limits p ^ and p ^ 1, only the 
edge channels (in Fig. 1(c), Figs. 3(a) and 3(c)) carry 
current and the point-contact conductance is bi-modal 
(see Table II). 

At criticality p = 0.5, the form of the point-contact 
conductance distribution P(Gpc) is more complicated. 
The distributions of the edge conductance obtained from 
numerical simulations of systems with Ly = 9 and var- 
ious Xd and Lx are shown in Fig. 6. Similar results are 
obtained for the distribution in the bulk (see Fig. 7). 
For Lx ^ Ly^ the distribution tends to a limiting form 
that depends on Xd, ^p, and Ly. For x^ ^ Ly, the 

dependence of this limiting distribution disappears 
(Fig. 7). Surprisingly, there is no self- averaging of the 
point-contact conductance even when x^ ^ 00 and the 
distribution remains broad. 

In Fig. 8, we show the squared fiux amplitudes 
|Qp(i = l,...,L°^*,...,R°^S...,iV). Note that in the 
asymmetric case, the current is distributed all across the 
sample even in the limit x^ ^ Ly (Fig. 8(a)). This is in 
sharp contrast to the symmetric case where the current 
quickly decays (Fig. 8(b)). 




(b) 













6 ^ 



0.2 



0.4 



0.6 



0.8 



Fig. 6. The distribution of the edge conductance for Ly = 9, (a) 
xa=S for Lx = 12(n), 36(o), 900(— ), and (b) x^ = 9 for L^ = 
24(+), 36(x), 900(— ). Ensemble averages over 1,000,000 systems 
have been taken. Irrespective of the value of x^, the dependence 
on Lx disappears for Lx ^ Ly. 
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Fig. 7. The distribution of the bulk conductance for Ly = 9, 
for = 9(0), 45(x), 450 ( — ). Ensemble averages over 5,000,000 
systems have been taken. Inset: Distribution of logarithms of 
point-contact conductance. For large x^, we see convergence to 
a broad distribution. 





Fig. 8. The squared flux amplitude |cip at criticality in (a) the 
asymmetric CC network of (Lx,Ly) = (900, 45) and (b) the sym- 
metric CC network of (Lx,Ly) = (900,44). Darker areas corre- 
spond to lower squared amplitudes. Probes, indicated by wedges, 
are attached at the middle of the system. Current flows from the 
left probe to the right probe (cj^ = 1, cJJ = 0, see eq. (10)). 



3.2 Dependence of (Gpc) on Xd 

To quantify how the conductance distribution con- 
verges to its hmiting form, we study the dependence 
of the averaged conductance. We have found that the 
averaged conductance converges exponentiahy, 



{Gpc) — {Goo) 



1 



6exp(-- 



(15) 



An example is shown in Fig. 9. Note that the values of 
{Goo)^ ^, and A depend, in principle, on and Ly. We 
emphasize that Goo = in the symmetric CC model. 




Fig. 9. The average of the bulk conductance as a function of 
for Ly = 9. The data are an average over an ensemble 
of 2,000,000 systems. The solid hue is the flt to eq. (15) with 
{G^> = 0.3208 ± 0.0001, b = 0.670 ± 0.017, and A = 4.67 ± 0.05 
(goodness of flt Q = 0.81). 



Scaling form describing the dependence of (Gpc) on 
Xd and Ly can be derived by assuming following factor- 
ization, 

(Gpc) = h{yp, Ly) f{xd, Ly). (16) 

To eliminate the ambiguity in this factorization, we set 
f{xd 00, Ly) = 1. Taking the limit Xd 00, 

(Goo) = %p,i,). (17) 

Comparing with eq. (15), we can write 

fixd^Ly) = l + 6exp(-Xd/A), (18) 

Xd=Xd/Ly, A=X/Ly. (19) 

We have found that data for different Xd and different Ly 
collapse onto a single curve (Fig. 10) with the following 
values, 



b = 0.675±0.019, A = 0.518±0.006. 



(20) 




3.3 Dependence of (Gpc) on Ly 

We now analyze the Xd-dependence of averaged edge 
conductance (Gpc) and similarly for the bulk conduc- 
tance (Gp^) for various Ly. 



Fig. 10. The ratio / = (Gpc) /(Coo) as a function of = x^jLy 
for Ly = 13(x), 25(0), 45(+) (edge conductance) and Ly = 
13(a), 25(n), 45(0) (bulk conductance). The sohd hue is a flt to 
eq. (18) (goodness of flt probabihty Q = 0.59). 
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Fig. 11. The imiting value of the point-contact conductance 
(Goo) for edge(x) and bulk(o) as a function of Ly. The data 
are estimated from disorder averages over 400,000 systems. Sohd 
Hnes are the fits to eq. (21). 



Table IIL Best-fit parameters for eq. (21). 





Co 


a 


Ci 


Q 


Edge 


0.6907±0.0053 


-0.3413±0.0016 


0.2531±0.0126 


0.59 


Bulk 


0.4590ib0.0018 


-0.2338ib0.0009 


0.4167ib0.0047 


0.70 



3.4 Dependence of (Goo) on Ly 

The dependence of (Goo) on Ly for edge and bulk is 
shown in Fig. 11. We have found that the foUowing form 



r^i T c 

^Q-Ly 



Ijy 



(21) 



fits our data. The best-fit values of parameters are listed 
in Table III. Here i denotes whether Gp^ (edge) or Gp^ 
(bulk). The first term is a non-trivial power law decay 
that reflects the multi-fractal^^"^ '^^ nature of the conduct- 
ing states. The second term is a correction for the dis- 
creteness of the model and the effect of the boundary. -^^^ 
The difference between edge and bulk conductance may 
originate from the difference between the surface and 
bulk multi-fractality.^''^ 

4. Summary and Concluding Remarks 

In this paper, we have calculated the point-contact 
conductance Gpc in the asymmetric Chalker-Coddington 
network model and found a novel metallic behavior. In 
contrast to the symmetric Chalker-Coddington network 
model, the point-contact conductance distribution con- 
verges to a broad distribution for a large separation of 
the contacts. This is true both for contacts attached to 
the bulk and the edges of the sample at criticality. This 
broad distribution reflects the nature of the current dis- 
tribution of the perfectly conducting state. We have also 
studied the averaged point-contact conductance in the 
limit of a large circumference, and found a scaling form. 
Both (G^) and (G^) show non-trivial power law decay 
eq. (21) (see also Fig. 11), which is a characteristic of 
criticality. 



So far we have focused on criticality p = 0.5. When 
p deviates from 0.5, the states are localized in the trans- 
verse direction. When system width Ly exceeds the 
transverse localization length, the perfectly conducting 
state is localized along one of the edges, 2/ = 1 for p > 0.5 
and y = Ly for p < 0.5 (Figs. 3(a) and 3(c) are ex- 
treme examples). In this case, (Gpc) decays quickly with 
the distance from the conducting edge. Broad distribu- 
tions P(Gpc) are observed only when we attach contacts 
near the conducting edge. Note that even in the ordinary 
quantum Hall effect, such fluctuations in point-contact 
conductances are expected near the edges. 

It is known that the conductance distribution is sen- 
sitive to the symmetry class (unitary, orthogonal, or sym- 
plectic) classified according to the presence or absence 
of time-reversal and spin-rotation symmetries. Since the 
time-reversal symmetry is broken in the scattering ma- 
trix eq. (7), the asymmetric Chalker-Coddington model 
belongs to the unitary class. A perfectly conducting 
channel also arises in the symplectic class, which is real- 
ized in carbon nanotubes.^^^^^ The distribution of the 
point-contact conductance in the symplectic class, espe- 
cially in the metal phase, may also be worth investigat- 
ing. 
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